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A primordial vector mode and its associated magnetic field generation are investigated. Firstly, 
we put an observational constraint on the amount of the primordial vector mode from the seven-year 
WMAP data. The constraint is found as r„ < 0.012, where r v and r are the amounts of 

O ' vector and tensor perturbation amplitudes with respect to the scalar one, respectively. Secondly, 

we calculate the spectrum of magnetic fields inevitably created from the primordial vector mode, 
. given the constraint on r v . It is found that the maximum amount of magnetic fields generated from 

the vector mode is given by B < 10 _22 G ( [) ' " 12 ) 1/ ' 2 ( 002 ) ^"" +1 ^ 2 with n v being a spectral index of 



the vector mode. We find a non-trivial cancellation of the magnetic field generation in the radiation 
dominated era, which creates a characteristic cut off in the magnetic field spectrum around k ~ 1.0 
Mpc" 1 . 

PACS numbers: 98.70.Vc, 95.30.-k, 98.80.Es 
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O ■ I. INTRODUCTION 

There are observational evidences which indicate that magnetic fields exist not only in galaxies, but also in even 
larger systems, such as cluster of galaxies and intra-cluster spaces [Tj. Yet, the origin of such large scale magnetic 
fields is still a mystery PJ . It is now widely believed that the magnetic fields at large scales are amplified from a tiny 
field and maintained by the hydro-magnetic processes, i.e., the dynamo. However, the dynamo needs a seed field to 
yjr^ ' act on and does not explain the origin of magnetic fields itself. As far as the magnetic fields in galaxies are concerned, 
the seed fields as large as 10 -20 ~ 10 -30 G are required in order to account for the observed fields Q of order l/i G 
C ■ at the present universe. Recent discovery of magnetic fields in galaxies at high redshifts [3] may require even larger 
seed fields. 

£SJ . A classical way to generate seed magnetic fields in astrophysics is based on the Biermann battery effect Q . The 
Biermann battery works in various astrophysical systems, such as stars supernova remnants M, Ml , protogalaxies 
Q, large-scale structure formation Q, and ionization fronts at cosmological recombination [Tol IllljT These studies 
show that magnetic fields with amplitude 10~ 16 ~ 10~ 21 G could be generated. However, the coherence-length of 
, seed fields generated by such astrophysical mechanisms may tend to be too small to account for galaxy-scale magnetic 
fields. 

On the other hand, cosmological mechanisms at the early inflationary epoch can produce magnetic fields with a 
large coherence length since accelerating expansion during inflation can stretch small-scale fields to scales that can 
exceed the causal horizon [I3 - [l5j . However, because no magnetic field is generated in simplest models with the usual 
electromagnetic field, it is necessary to introduce some extensions to the standard particle model. Furthermore, it 
is recently argued that the backreactions from the electro-magnetic fields will stop the inflation and significantly 
suppress the magnetic field generation [IiKIt]], if they are properly taken into account. 

Generation mechanisms of magnetic fields in a decelerating universe prior to cosmological recombination have also 
been proposed. Originally, Harrison found that the vorticity in a primordial plasma can generate magnetic fields [l8j . 
This is because electrons and ions would tend to spin at different rates as the universe expands due to the radiation 
drag on electrons, arising a rotation-type electric current and thus inducing magnetic fields. Following his idea, many 
authors have investigated magnetic field generation through the second order vorticity generated from the first order 
density perturbations [l9l - l27| . 
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In this paper, we consider a first order vector (vorticity) mode firstly investigated by Rebhan followed by Lewis 
[1^,[2^]. In a Friedmann universe with a perfect fluid without anisotropic stress, the vector mode has only a decaying 
mode. This means a diverging vector metric perturbation (frame-dragging potential) at the initial time and therefore 
the model is inconsistent with an almost isotropic Friedmann universe, and goes beyond the linear perturbation 
theory. However, in the existence of anisotropic stress by free streaming particles such as neutrinos, it has been 
found that there exists a regular (growing) mode with an initial non-zero vorticity and with isotropic initial phase 
space distributions. We first present an observational constraint on this primordial vector mode amplitude using the 
seven-year WMAP data. We then estimate the amplitude of the magnetic fields inevitably created from this vector 
mode in the light of the obtained constraint on the vector mode amplitude. 

This paper is organized as follows. In Sec. II, the overview of the primordial vector mode investigated by [2^, [29[ 
is described and we give an analytic solution obtained by a tight coupling approximation. In Sec. Ill, we show the 
constraints on the vector mode amplitude and its spectral index from the seven-year WMAP data and discuss its 
implication. In Sec. IV, we calculate the magnetic field spectrum inevitably generated from the primordial vector 
mode. Finally, Sec. V is devoted to the summary and discussion. 



II. EQUATIONS AND SOLUTIONS 



A. preliminary 

Here we quickly review basic equations for the evolution of primordial vector modes. We consider linear perturba- 
tions in the synchronous gauge, in a flat Friedmann-Robertson- Walker universe with a metric 

ds 2 = a{rj) 2 [-drj 2 + (5 y + h^dx 4, dx>] , (1) 

where a(j]) is a scale factor. Because we are interested in the vector mode, the metric perturbations are decomposed 
in Fourier space as 

= ikihj + ikjhY , (2) 

where bX is divergence-less, i.e., k l bX = 0. It is useful to expand the vector bX in terms of two independent transverse 
basis vectors as 

hY(k) = ^2hl(k)et(k), (3) 
± 

where ± represents the parity. Because each parity component and each Fourier mode evolve independently, we omit 
k and ± dependencies for simplicity in the rest of the paper. Let us work with the gauge invariant metric variable, a, 
which is defined using the synchronous gauge variable h as a = hX jk. Then, the linearized Einstein equations give 

k 2 a = -IGnGa 2 (p + P) v , (4) 

a + 2U(j = &TrGa 2 PTr/k , (5) 

where dot denotes a derivative with respect to the conformal time r/, H = a/a is the conformal Hubble parameter, 
p and P are the zero-th order density and pressure of the total fluid, respectively, and v and 7r are the velocity and 
anisotropic stress in the vector mode, respectively. From Eq.© it is readily found that the vector metric perturbation 
has only a decaying mode a oc a~ 2 if there is no anisotropic stress. 
For the matter part, the Euler equation of baryon is given by 



ib + T-Lv b = --^-an e OT (vb - u 7 ) , (6) 
opb 



where n e is the electron number density and cry is the Thomson scattering cross section. For the photons we expand 
the distribution function for the vector mode into multipoles and rewrite the Boltzmann equation as 

Uy + \kit 1 = —an e aT (f 7 — v b ) , (7) 
8 

8,8, / 9 9\8, 

7r 7 + ~kh - -kv 7 = -an e a T I — tt 7 - ~E 2 I - -ka , (8) 

where In is the £-th order angular moment of the photon distribution, and Ei is the £-th order moment of the E-modc 
polarization (29[. 
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B. tight coupling approximation and solutions 



In the very early universe photons and baryons are tightly coupled as the opacity A = an e aT is large. This enables 
us to expand the equations by a tight coupling parameter, 



T 



i 



0.02 



(10) 



where is the baryon density normalized by the critical density, and h is the Hubble parameter normalized by 100 
km/s/Mpc. At the lowest order, the equation for photon velocity is 
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where R = f^. The solution is 
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where u 7 .i n i is the initial photon velocity. Therefore, the fluid velocity stays constant deep in the radiation dominated 
era where fic 1. At this order the baryon-photon slip and the photon anisotropic stress do not exist. At next order 
they are given by 
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respectively, where is the solution of the vector metric perturbation at the lowest order. Note that the transfer 
of baryon-photon slip is independent of k. In other words, because R rj R/t the slip term in the vector mode is 
further suppressed by a factor of Rj (fcr) on small scales. This is compared with the slip term in the scalar mode, i.e. 
w 7 — Vb oc (k/f)R, where the k dependence derives from the pressure gradient. 

For the later discussion about magnetic field generation we derive the baryon-photon slip at the next order in the 
tight coupling approximation. We find it to be 
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Therefore, up to the second order in the tight coupling approximation the baryon-photon slip is given by 
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where we have neglected cosmological redshift terms. Note that the signs of the two terms are different. We shall 
see below that, because of this, a significant amount of magnetic fields resolves away when the two terms become 
comparable around the Silk damping epoch. 

Before moving to the next section we summarize the initial conditions [28l [29J. By expanding the equations in 
powers of kr and assuming the radiation dominated era, we find at the lowest order 
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where we have defined the ratios R v = — Q 

Pu-t 

in Fig. H 



and R-y = 1 — R v . Evolutions of the vector perturbations are presented 
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FIG. 1: Time evolutions of the vector perturbation variables with wavenumber k — 0.1 Mpc -1 (left) and k = 1.0 Mpc -1 (right). 
The vector potential a starts to decay after horizon crossing, while velocities Vb,-y stay constant until the Silk damping takes 
place. 



III. CMB CONSTRAINT ON THE PRIMORDIAL VECTOR MODE 



In order to compare with the CMB data precisely, we use the CosmoMC package with a modification to include 
two new parameters, r v and n v . Here r v is the vector-scalar ratio and n v is the spectral index of the vector power 
spectrum. Specifically, we parameterize the power spectrum of primordial vector mode with an initial amplitude of 
the metric perturbation a as 

^^a(k)a%k'))=V a (k)8(k-k')=AJ^-\ 5(k-k'), (22) 

where ko = 0.002 Mpc -1 is the pivot scale. The vector-scalar ratio is defined by 




(23) 



where A s is the scalar counterpart of the power spectrum amplitude. 

The likelihood function we calculate is given by the WMAP collaboration [3(J|3l| and the codes are publicly available 
at their web site. To include the primordial vector mode, we calculate the angular power spectrum, C ve ctor i, of the 
CMB anisotropy of temperature and E-mode polarization auto-correlations (TT and EE) and their cross correlation 
(TE) using the CAMB code [32]. We add these power spectra to those from the scalar and tensor modes as 

^yrr,TE,EE _ .^tt,te,ee . ^tt, te.ee . ^.tt, te.ee , a \ 

°tot I ~ ° scalar I ' °tonsor t ' ° vector t > \ LqL ) 

where we have assumed that the primordial scalar, vector and tensor modes are statistically independent. Then Ct t I 
are fitted to the CMB data. The B-mode polarization angular power spectrum is not used for our analysis because 
the sensitivity of the current observational data of B-mode polarization is not enough to give limits on the vector or 
tensor mode amplitudes. 

In Fig [2l we depict the CMB angular power spectra along with the seven-year WMAP data. In that figure, the 
parameters of the primordial vector mode are fixed to r v = 9.45 x 10 -3 and n v — 0.921, respectively, which are the 
allowed values at 95% confidence levels (red dashed lines in Fig|2]) . The spectra look similar to those from the tensor 
mode perturbations, although the powers extend to the higher multipolcs without oscillatory features. As clearly 
seen from the figure, the current constraints are mainly placed on from the TT power spectrum at low multipoles if 
the vector mode is nearly scale invariant. The situation is the same with the case of the current constraint on the 
primordial tensor mode (gravitational waves) from the CMB power spectrum. Only in the panel for the TT power 
spectrum we also show the case with a bluer spectral index n v = 2.0 with r v — 1.0 X 10~ 3 , which are also at the edge 
of 95% confidence levels. In this case, we find that the constraint again comes from TT power spectrum, but at the 
higher multipoles around £ « 1000. 

The correlation between the parameters of primordial vector modes, r v and n v , is shown in the right panel of Fig 
[3l If the primordial vector perturbation is given almost scale invariant, we found that the constraint is put at lowest 
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FIG. 2: CMB angular anisotropies power spectra, temperature (TT), E-mode polarization (EE), and their cross-correlation 
(TE), induced by scalar and vector type perturbations. The constraint mainly comes from the TT angular power spectrum. 
The vector mode parameters are taken as (r„,n„) = (0.945 x 10~ 2 , 0.921) and (1.0 x 10~ 3 ,2.0) for red dashed lines and green 
dash-dotted line (TT spectrum only), respectively. 



multipoles as mentioned above and the constraint is the weakest. When the spectrum becomes bluer as n v > 1, the 
constraint comes from the higher multipoles and hence it becomes tighter. For example, when n v « 2.0 the constraint 
becomes r v < 0.001, while r v < 0.009 when n v fts 1.0. 

To make our analysis as general as possible, we made our Markov chain analysis with and without the primordial 
tensor mode. The effect of including the tensor mode on the constraint on the vector mode is also shown in the right 
panel of Fig. [3] We observe that the constraint on the vector mode generally becomes tighter if the tensor mode is 
included. The constraints we found are r v < 9.45 x 10~ 3 (without tensor) and r v < 8.36 x 10 -3 (with tensor) when 
the all the other parameters are marginalized. The simultaneous constraint on the vector and tensor modes is shown 
in the left panel of Fig. [3] The result shown is easily understood: because the angular power spectra of temperature 
anisotropies from the vector and tensor modes are similar, the current CMB data can put constraint on the total 
amount of the vector and tensor mode perturbations. We found the constraint to be 

rv + (^) < 0.012 , (25) 

at 95.4% confidence level where r stands for the tensor-scalar ratio. 



IV. MAGNETIC FIELD GENERATION 



When there exists velocity difference between baryons and photons in the vector mode, magnetic fields may be 
generated inevitably. This arises because photons scatter off electrons preferable to protons. Electric fields are 
induced to prevent charge separation between electrons and protons, and magnetic fields are generated from these 
induced electric fields through Maxwell equations if there exists the vorticity difference. In general, magnetic fields 
then affect the evolution of the fluids. However, we omit any backreactions from magnetic fields on the fluid motion 
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FIG. 3: 68 % and 95 % confidence regions from WMAP 7-yr data on r-r v (left) and n v -r v planes (right). Left panel shows 
that constraint is on the total amount of vector and scalar components. Right panel shows constraints become severer when 
vector perturbation spectrum gets away from scale invariant n v — 1. In that panel red contours are obtained from the analysis 
without tensor modes, blue contours are obtained including tensors. Blue and red spectra are allowed when r v —¥ 0. 



because we are interested in the magnetic field generation from zero and thus the backreactions from the magnetic 
fields will be negligible. In terms of the cosmological perturbation theory, the backreactions from magnetic fields are 
second order, namely, the Lorentz force would be vB ~ 0(v 2 ) and the energy momentum tensor will be proportional 
to B 2 ~ 0{v 2 ). 

The evolution equation of magnetic fields due to the Thomson scattering is given by 

d(a 2 B l ) 4cr T p 7 a ilk 

dt = ^f- £ K A* " Vb o,k) , (26) 

where e y& is the Levi-Civita tensor and e is the electric charge. For the scalar type perturbation, this term vanishes 
because e l ^ k v^ k ~ e l ^ k kjkkV S — 0. In this paper we solve the above equation with the initial condition B = at 
z = 10 9 , which roughly corresponds the time of neutrino decoupling. 

Even in the case where the vector mode perturbation exists, the mean of the magnetic fields should still be zero, 
but there exists the variance. To find the spectrum of magnetic fields, we first calculate B l Bi. That is given by 

a 4 B l (k,t)B*{k>,t) = {8 jl 8 km -8 jm 8 kl )k k k' m J J 'o(0*«i(j?,< / )o(*")*«/(^0*'*" ■ ( 27 ) 

Taking the ensemble average of this expression and defining the vector power spectrum as 

(Sv 3 (k,t')8v* e (k>,t")^ = ^-V a (k)P jt (k)Sv(k,t')Sv{k,t")S{k - k') (28) 

P je (k) = Sji-kjk , (29) 

where 8v(k,t) is the transfer function of baryon-photon slip with a = 1 at the initial time, and the magnetic field 
power spectrum can be written as 



a 



4 (t)^S B (k,t) = (^-) 2P t7 (k)k 2 



2 



dt'a(t')p-y(t')Sv(k,t') 



2 



(30) 



2n 2 ° y ' ' \ 3e J 
Here, magnetic field power spectrum Ss{k) is defined as 

S B {k)8{k-k') = ^B l (k)B z {k')) . (31) 

We numerically calculate the magnetic field spectra and show them in Fig. |U In that figure, we took the vector mode 
parameters as r v = 0.01 and n v — 1, which give the maximum amount of the vector mode allowed at large scales. 
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First of all, let us roughly estimate the amplitude of the generated magnetic fields. From Eq. ([26|) . it is estimated as 



(« B ) - ( -3^- j JjPiK v i - v b) i 

/4a T \ 2 i? /4i?„ + 5\ 



Vaer^UJ" (32) 

where H is the usual Hubble parameter and we have used the tight coupling solution for (u 7 — Vb) derived in Sec. II. 
By substituting Eq.flTD), p 7 ~ 2 x HT 51 (1 + z) 4 [GeV 4 ], and a T ~ 1.7 x 10 3 [GeV" 2 ], we find 

.' B ~1.2X10-»GC * U^V'f^V- (33) 



VMpc^V V !0 4 / V0.01 

Therefore, around the cosmological recombination epoch, B ~ 10~ 21 G is expected for the magnetic field strength. 

At super-horizon scales, the magnetic fields have a power-law B oc k if the vector mode power spectrum is scale 
invariant. This is manifestly shown by Eq. (|32p . If the primordial vector mode is tilted as n v ^ 1, the magnetic field 
spectrum at super horizon scale should become as S oc k^""" +1 '' 2 . 

At subhorizon scales, we found some interesting features in the resultant magnetic field spectrum. First, in the 
radiation dominated era, there is a peak just above the Silk damping scale. Secondly, we find a characteristic cut-off 
below the Silk damping scale. These features were not observed in the magnetic field spectrum from the second 
order vector modes. In that case the magnetic field spectrum is extended toward much smaller scales [24H26j |. In the 
primordial vector mode considered here there is no significant difference in the evolution of velocity difference between 
photons and baryons before and after the horizon crossing. Therefore the magnetic fields continue to grow after the 
horizon crossing, which is manifestly shown in the magnetic field spectrum as B oc k for k > kh, with kh being the 
wavenumber corresponding to the Hubble scale. It is expected that the generation of magnetic fields ceases around 
the Silk damping epoch where the velocity difference, i> 7 — Vb, reaches its maximum value and starts to diminish. The 
peak position and the amplitude of magnetic fields there can be estimated as follows. The magnetic fields generated 
from Ea.(|26[) is estimated as 

B = {a 2 B) r; — p 7 /c(w 7 — Vb) oc ka , (34) 
H 

where we have defined the comoving magnetic field B, and used the relations in the radiation dominated era such 
as H oc a~ 2 , p 7 oc a~ 4 , and (w 7 — Vb) oc a 2 . The comoving magnetic fields evolve in time as oc a and therefore the 
maximum value is reached at the Silk scale where the perturbations start to be erased. Because the diffusion scale is 
scaled as kdis oc a -3 / 2 [33[, the scale factor when diffusion damping occurs at a given scale k is given as asuk oc fc~ 2 / 3 . 
Therefore, if the magnetic field generation just ceased when the Silk damping started we should expect that the 
magnetic fields have the spectrum as B oc kasm oc A; 1 / 3 . In Fig. 0] we find that the peak amplitudes at different 
redshifts are indeed on this scaling relation. 

However, we found a non-trivial cancellation of the magnetic field generation. In fact, the magnetic field generation 
does not only cease at the Silk damping epoch, but a significant amount of magnetic fields generated by that time 
vanishes as is shown in the right panel of Fig. 5. This happens because the sign of the velocity difference, i> 7 — i>b, flips 
at the onset of the Silk damping, which is shown in the left panel of Fig. 5. Deep in the radiation dominated era, the 
velocity of radiation (photon) fluid stays constant due to the redshift of energy density (Eq. (fl2"|) ). while the velocity 
of baryon fluid should decay as x>b oc a -1 if baryons had no interaction with photons. Therefore during that time the 
baryon fluid is dragged by the photon fluid, and we expect that v 7 — w& > if w 7i i n i > 0. After the diffusion damping 
starts to erase the perturbations in photons, the baryon fluid, which will try to keep rotating by its inertia, drags the 
photons so that w 7 — Vb < if w 7i i n i > 0. This argument is consistent with the result obtained from the tight coupling 
expansions, in other words, the magnetic fields generated by the first term in Eq. (1161) are largely compensated by 
the second term. The resultant (comoving) magnetic fields spectrum has the characteristic cut-off at k > 1 [Mpc" 1 ]. 



V. SUMMARY AND DISCUSSION 

In this paper we examined observational constraints on the primordial vector mode discussed in [29| . This vector 
mode is a vector analogue of neutrino isocurvature velocity scalar mode ,34| and has a non-decaying solution of 
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fc[l/Mpc] fc[l/Mpc] 

FIG. 4: Left: The spectra of magnetic fields generated from primordial vector modes at several times as indicated in the figure. 
The final spectrum at decoupling epoch is shown by the black solid line. Right: The dependence on n v at z = 1 100. 
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FIG. 5: Time evolutions of the source for the magnetic field (left) and the transfer of the magnetic field (right) for some 



interactions between photons and baryons, and therefore photons drag the baryons through Thomson interactions (the first 
term in Eq. JT6J)). 



vector metric perturbation. The vector mode perturbations generate CMB angular anisotropies mainly through the 
Dopplcr effect, whose angular power spectrum is similar to that from the tensor mode without oscillatory features 
at small angular scales. Thus the amplitude of vector mode perturbations is constrained from CMB data. We used 
seven-year WMAP data to place a constraint on the vector-scalar ratio r v . The constraint we obtained is r v < 0.01 
when the vector spectrum is nearly scale invariant. If the tensor mode is included the constraint is generalized to 
r v + (r/40) < 0.012 where r is tensor-scalar ratio. The current constraint is placed mainly from the temperature 
anisotropy data. Future precise polarization data will tighten this bound significantly, down to r v < 10~ 3 for Planck 

M 

As discussed by several authors, magnetic fields are generated before cosmological recombination through 
Compton scatterings if there existed cosmological vector mode perturbations. We numerically calculate the 
magnetic field spectrum generated from such primordial vector mode. We found that the spectrum is B rs 

10~ 23 G (o%-)^ 2 ((TuTs) " " +1 ^ 2 at cosmological recombination where n v is the spectral index of the primordial vector 
mode. The magnetic fields monotonically increases with wavenumber as B oc fcC"^ 1 )/ 2 U p to the Silk damping scale. 
Below this scale we found that there is a non-trivial cancellation in the magnetic field generation, leading around five 
orders of magnitude decrease in the magnetic field amplitude. The cancellation occurs because the baryons drag the 
photons around the Silk damping epoch, before which, on the contrary, the photons dragged the baryons. 

In fact, perhaps surprisingly, using the second order tight coupling approximation we can show that the magnetic 
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fields should eventually vanish at this order as follows. Deep in the radiation dominated era and at sub-horizon scales, 
the tight coupling equation of photon velocity is written as 



k 2 4 

«; + T -« 7 = 0. (35) 



This equation can be solved to give the solution 



v 1 =v ltini exp I -— 77 I , (36) 

where we have introduced a constant /3 to describe the differential optical depth as f = /3r)~ 2 in the radiation 
dominated era. Inserting this solution into the tight coupling expression of the baryon-photon slip fEq. ([TE| ). and 
writing the time dependences of a, R and p 1 as a — ar), R = R^rj and p 1 — JJ]^ 4 ", respectively, the time integral to 
get the magnetic field amplitude can be expressed as 



4 12a 3 ( 4fc2 3 



dr] . (37) 



The above integral of the form J^ md (1 — Krf 3 ) exp(—Kr] 3 /3)dr] = [77 exp (— -| if 77 s )] ^ cnd (with K constant) goes exactly 
to zero for ?7i n i — > and 77 en d — > 00. 

The cancellation is perfect only within the second order tight coupling approximation in the completely radiation 
dominated era, and only when we integrate from r\ = to 77 = 00. Because we assume that the magnetic field 
generation starts at some finite time, B(rji n i) = with ?7i n i = i](z = 10 9 ), the surface term should remain as 

a B ~ - I 3g I a 7 I — I v 7 , in i?7im = - I 3e I a /9 7 yV 7 , in i • (38) 

We have found that the surface term above dominates the magnetic fields at small scales over the higher order terms 
in the tight coupling approximation, and indeed, our numerical solution is consistent with the above estimate at 
k > 10 Mpc -1 . This suggests that the magnetic fields start to increase again with wavenumber as B oc k^ n,J+1 ^ 2 
which is shown in Fig. 01 Consequently, magnetic fields will have a small-scale power up to ~ 2 x 10 -3 pc in the 
comoving scale. This scale corresponds to the damping scale when the temperature of the universe was around MeV, 
after which epoch neutrinos can free-stream and the primordial vector mode solution considered here can apply. 

The vector modes and magnetic fields considered in this paper are different from the recently investigated second 
order vector modes or magnetic fields (3^413 • In those studies the authors considered the second order vector modes 
and/or magnetic fields generated from the non-linear couplings of the first order scalar (density) modes. In the 
framework of the cosmological perturbation theory, the second order solution may be considered as a particular 
solution (in the limit of neglecting vector-scalar couplings), because the couplings of first order density modes can 
be considered as a source term for the otherwise homogeneous evolution equations of the vector mode. Therefore, 
the general solution may be expressed as a superposition of the solutions, namely, the second order solution and the 
solution considered in this paper. The second order solution suggests that magnetic fields generated in the second 
order vector modes have the amplitude asB~ 10 -26 G around recombination at k k, O.OlMpc^ 1 . Comparing the 
amplitudes, the magnetic fields from the homogeneous solution considered in this paper will dominate the density 
perturbation induced magnetic fields if r v > 1CP 8 . 

Another way to regularize the cosmological vector mode in the early universe will be a modification of the vector 
sector of gravity. An interestin g^ examp le is Einstein- Aether theory which has been recently investigated with cosmo- 
logical perturbations (see, e.g., [38|, |39(). In this class of models the vector metric now becomes a dynamical variable, 
and the Aether field can induce another growing mode in the vector cosmological perturbations. It is expected that 
a sizable amount of vector type velocities in the fluids is induced [39| . The magnetic field generation from this vector 
mode will be an interesting subject, but beyond the scope of this paper. 

In conclusion, in this paper it is found that the magnetic fields generated from the primordial vector mode at 
recombination are allowed at most B < 10 -21 G at k « 0.1 Mpc -1 . The upper bound comes from the upper bound 
on the primordial vector mode amplitude obtained from the seven-year WMAP data. This field amplitude is too small 
to be directly observed through Faraday effects on the distant radio sources and/or delays of high energy photons 
from gamma-ray bursts or blazars, however it may be large enough for the fields to be a seed for galactic magnetic 
fields observed today. 
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